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Abstract. The hyperdeterminant of format 2x2x2x2 is a polynomial 
of degree 24 in 16 unknowns which has 2894276 terms. We compute the 
Newton polytope of this polynomial and the secondary polytope of the 4- 
cube. The 87959448 regular triangulations of the 4-cube are classified into 
25448 Z?-equivalence classes, one for each vertex of the Newton polytope. 
The 4-cube has 80876 coarsest regular subdivisions, one for each facet of the 
secondary polytope, but only 268 of them come from the hyperdeterminant. 



1. Introduction 

This article makes three contributions to the mathematics of computation: 

• solution to a computational challenge problem posed by I.M. Gel'fand, 

• new theorems in algebra (hyperdeterminant) and geometry (4-cube), 

• new computational methodology (algorithms and software) for large 
polynomials and large convex polytopes in the presence of symmetry. 

We begin our discussion with a small familiar polynomial in four unknowns: 

D22 = CooCii - CoiCio- 

This irreducible quadratic polynomial is the determinant of a 2 x 2- matrix (cj-, ). 
Next, the hyperdeterminant of a 2 x 2 x 2-tensor is the following irreducible 
polynomial in eight unknowns which is the sum of twelve terms of degree four: 

D222 = 4coooCoilCioiCiio + 4cooiCoioCiooClll 

+Cooo'^lll + CpoiCfj^Q + CqiqCioi + CoiiCfoo 
— 2CoooCooiClloClll — 2CoooCoioClOlClll — 2CoooCoilCiooClll 
— 2cooiCoioCioiCiio — 2cooiCoilCiioCioo — 2coioCoilCioiCioo- 

This polynomial is known as the tangle in the physics literature. Hyperdetermi- 
nants of arbitrary format were studied by Gel'fand, Kapranov, and Zelevinsky 
|12j . They have numerous applications ranging from quantum information 
theory [T31[T3] to computational biology [HHI] and numerical analysis [7| 1^. 
Basic definitions concerning hyperdeterminants will be reviewed in Section 3. 

The starting point of this paper is a problem stated by I. M. Gel'fand in his 
Fall 2005 research seminar at Rutgers University: find the monomial expansion 
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of the hyperdeterminant -D2222 ofa2x2x2x 2-tensor. Back in 1852, Schlafli 
had given a nested formula (see §14.4]) for D22221 which is a polynomial 
of degree 24 in 16 unknowns Coooo; Cqooi; • • • ; Cim, and, 150 years later, Luque 
and Thibon ^B] expressed -D2222 in terms of the fundamental tensor invariants. 
Gel'fand's challenge to the computer algebra community was to expand these 
formulas into monomials. We solved Gel'fand's problem: 

Theorem 1. The hyperdeterminant -D2222 "is the sum 0/ 2894276 monomials 
in 9617 orbits. The Newton polytope of D2222 'is an 11- dimensional polytope 
with 268 facets in 8 orbits, and 25448 vertices in 111 orbits. It contains the 
exponent vectors 0/ 20992 monomials in 69 orbits which do not appear in -D2222- 

In this theorem, the term orbits refers to the full symmetry group -B4 of 
the 4-dimensional cube, which has order 384, and the Newton polytope is the 
convex hull in M}^ of the exponent vectors of all monomials appearing in -D2222- 

Our computational proof of Theorem 1 is explained in Sections 3 and 4. We 
list representatives for each orbit of facets, vertices, and "missing monomials," 
and we discuss key properties found in these data. Complete data are available 
at our supplementary materials website bio.math.berkeley.edu/4cube/ 

According to the general theory of Gel'fand, Kapranov and Zelevinsky ^21 
§11], there is a beautiful and deep relationship between the monomials in the 
2 X 2 X ■ ■ ■ X 2-hyperdeterminant and the triangulations of the n-dimensional 
cube. In the language of polyhedral geometry, this relationship can be stated 
succinctly as follows: the secondary polytope of the n-cube equals the Minkowski 
sum of the Newton polytopes of the hyperdeterminants corresponding to all 
the faces of the n-cube, up to dimension n. Readers unfamiliar with these 
geometric concepts will find an elementary and self-contained introduction in 
Section 2, where we explain the correspondence between the twelve monomials 
in D222 and the 74 triangulations of the 3-cube. 

Theorem 1 opened up the possibility of determining the same correspondence 
for the 4-cube. However, the challenges were formidable, in light of the size of 
the secondary polytope of the 4-cube: 

Theorem 2. The A-cube has 87959448 regular triangulations in 235277 orbits. 
The secondary polytope of the A-cube has 80876 facets in 334 orbits. 

The first statement in Theorem 2 is a known result due to Pfeifle and Ram- 
bau [inj, who had computed all regular triangulations of the 4-cube using 
Rambau's software TOPCOM We independently verified their enumeration. 
The second statement in Theorem 2 is one of the contributions of this paper. 

Equipped with Theorems 1 and 2, it took us two more months of serious 
computational work to obtain what we consider to be the main mathematical 
result of this paper: the extension of the material in Section 2 from dimension 
three to dimension four. Each regular triangulation of the 4-cube maps to a 
monomial in -D2222, and the fibers of this map are the D -equivalence classes 
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of [12]. We found that, up to symmetry, there are 111 /^-equivalence classes. 
They are discussed in Section 6. Section 7 studies the 80876 facets in Theorem 
2, and it classifies the coarsest regular polyhedral subdivisions of the 4-cube. 
Most of the algorithms and softwares we developed and used here can easily be 
generalized and applied to other polynomials and polytopes with symmetry. 

Section 5 places our computations into a larger mathematical context. The 
hyperdeterminant -D2222 is the A- discriminant for the matrix of the 4-cube 



A 



1111111111111111 
0000000011111111 
0000111100001111 
0011001100110011 
0101010101010101 



Dickenstein, Feichtner and Sturmfels [8^ recently gave a recipe for computing 
the Newton polytope of the A-discriminant for any point configuration A. Our 
computations are complementary to their tropical approach. We represent 
subdivisions of the cube A by their dual polyhedral complexes, here called 
tight spans as in The related cases when ^4 is a second hypersimplex or 
a product of two simplices correspond to finite metrics j20j and to tropical 
convexity P respectively. We also study the irreducible factorizations of all 
leading forms of -D2222 that are supported on facets of the Newton polytope. 



2. The secondary polytope of the 3-cube 

As a warm-up for our study of the 4-cube, we first discuss our primary 
objects of interest for the 3-cube. Recall that the Newton polytope A/'(G) of 
a polynomial G is the convex hull of the exponent vectors of the monomials 
which appear in the expansion of G. The Newton polytope M{D222) of the 
hyperdeterminant -D222 is the convex hull in of the six rows of the matrix 
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These six monomials labeling the rows are the extreme monomials of -D222; 
which means that their exponent vectors are vertices of the Newton poly- 
tope M{D222)- The other six monomials in -D222 are not extreme monomials, 
since their exponent vectors lie in the relative interiors of edges of the four- 
dimensional polytope M{D222)- Combinatorially, ^[{0222) is a bipyramid over 
a tetrahedron. Its f -vector records the number of faces of dimensions 0, 1, 2, 
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and 3 respectively: 

(1) f{^f{D222)) = (6,14,16 , 8). 

The Newton polytope of the 2 x 2 x 2-hyperdeterminant has the following 
irredundant presentation by linear equations and inequalities: 

M{D222) = {(xooo,a^ooi, • • • e M*^ : 

2^000; ^ooij ^010^ 2;oii, Xioo, a^ioi, Xuq, Xiu > 0, 

3^000 + ^001 + Xqio + Xoil = Xqqq + Xqoi + a^lOO + ^101 = 2, 
^000 + ^010 + ^100 + ^110 = ^001 + ^011 + ^101 + ^111 = 2 }. 

In Section 4, we determine the analogous presentation for N {02222)- 

We consider the principal determinant oi the 3-cube. By jT2| Thm. 10. B. 1.2], 
this is the following product of determinants associated to faces and vertices: 

E222 = D222 ■ (coooCoii — cooiCoio) ■ (coooCioi — cooiCioo) ■ 

(coooCllO — CoioCloo) ■ (cooiClll — CoilCioi) • 
(CoioClll — CoiiCiio) ■ (ClooCiii — CloiCiio) • 

Cooo ■ cool ■ colo ■ Coil ■ Cloo ■ Cloi ■ Clio ■ Cm. 

The expansion of this polynomial of degree 24 has 231 monomials, of which 
74 are extreme monomials. The symmetry group of the 3-cube, which is the 
Weyl group of order 48, acts on these 231 monomials. The 74 extreme 
monomials come in six orbits: 

Type Monomial GKZ Vector Orbit Size 

1 -4coooc^oiC^ioCoiicfooCioiCiiocfii (15515115) 2 

2 -coooCooiCoioc|]iic5ooCioiCiiocfii (14436114) 8 

3 coooCooiCoioCoiiCiooC?oiCiioCiii (13446 213) 24 

4 CoooCooiCoioCoiiCiooCioiCiioCiii (133 5 5331) 12 

5 ~CoooCqqj^Cq;^qCq;^j^C]^qqC]^q]^C]^;^qC]^;^]^ (133 56222) 24 

6 '^ooo'^ooi''Oio''Oii''ioo'^ioi'^iio'^iii (2 226622 2) 4 

The Newton polytope M{E222) of the principal determinant E222 is the sec- 
ondary polytope of the 3-cube. It is 4-dimensional and has the /-vector (74, 152, 
100, 22). A Schlegel diagram of this polytope was posted by Julian Pfeifie at 

www . eg-models . de/models/Discrete_Mathematics/ 
Polytopes/Secondary_Polytopes/2000 . 09 . 031/ . 

Using POLYMAKE [TT], we found that the secondary polytope of the 3-cube 
has the following irredundant presentation by linear equations and inequalities: 



THE HYPERDETERMINANT AND TRIANGULATIONS OF THE 4-CUBE 



5 



Af{E222) = {(xooo,2;ooi, ••• ,3^111 ) e : 

3^000 + 3^001 + ^010 + ^011 = ^000 + ^^OOI + ^^lOO + ^^lOl = 12, 

Xooo + ^010 + ^100 + ^110 = ^001 + ^011 + ^101 + ^111 = 12, 
1 < Xqqo, Xqqi, Xqio, a^oii) 2;ioo, a^ioi, Xuq, Xm < 6, 
3^000 + ^001 ^ 4, xqoo + 3;oio ^ 4, xqoo + ^100 ^ 4, 
a^ooi + xou > 4, xoio + xqu > 4, xqoi + xwi > 4 }• 

According to the theory of Gel'fand, Kapranov, and Zelevinsky ^2); the 74 
vertices of the secondary polytope A/'(-E222) are in bjiection with the regular 
triangulations of the 3-cube. For the 3-cube all triangulations are regular, 
so the 3-cube has 74 triangulations which come in 6 orbits. The vertex of 
M[E222) corresponding to a given triangulation 11 is called the GKZ vector of 
n. The coordinate of the GKZ vector of 11 is the sum of the normalized 
volumes of all tetrahedra in 11 which contain the ith vertex of the cube. For 
a general introduction to triangulations, their GKZ vectors, and secondary 
polytopes, see the forthcoming book by De Loera, Rambau and Santos [TUj . 
The six types of triangulations of the 3-cube are depicted in [TUl Figure 1.38] 
and Figure 3]. Their tight spans, which are polyhedral complexes dual to 
the triangulations, are shown in Figure 1 below. Here is a detailed description: 

Type 1: These two triangulations divide the cube into five tetrahedra. There 
is one central tetrahedron of normalized volume two, and four unimodular 
tetrahedra which border the central tetrahedron. 

The remaining five types of triangulations divide the 3-cube into six uni- 
modular tetrahedra. Each of these triangulations uses a diagonal of the cube. 
Type 2: We get these 8 triangulations by shcing off the three vertices adjacent 
to a fixed vertex. The remaining bipyramid is cut into three tetrahedra. 
Type 3: These 24 triangulations also use a diagonal. Of the other six vertices 
we pick two that are diagonal on a facet, and we slice them off. 
Type 4: These twelve triangulations are indexed by ordered pairs of diagonals. 
The end points of the first diagonal are sliced off, and the remaining octahedron 
is triangulated using the second diagonal. 

Type 5: These 24 triangulations are indexed by a diagonal and one other 
vertex. That vertex is sliced off, and the remaining polytope is divided into a 
pentagonal ring of tetrahedra around the diagonal. 

Type 6: These four triangulations are indexed by the diagonals. The cube is 
divided into a hexagonal ring of tetrahedra around the diagonal. 

We depict each triangulation of the 3-cube by its corresponding tight span, 
which is the planar graph dual to the triangulation. Each vertex of the tight 
span corresponds to a tetrahedron in the triangulation, and two vertices are 
connected by an edge if the corresponding tetrahedra share a triangle. Regions 
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Figure 1. The tight spans dual to the six types of triangu- 
lations of the 3-cube. Each tight span represents an orbit of 
extreme monomials of the principal determinant -£'222- 



of the tight span correspond to interior edges of the triangulation. The six 
types of tight spans are shown in Figure 1. Compare it to Figure 2 in fH^ . 

The 22 facets of A/'(-E222) correspond to proper subdivisions of the 3-cube 
which are as coarse as possible. The eight facet inequalities like xqoo > 1 
correspond to slicing off one vertex. The eight inequalities like Xqoo ^ 6 cor- 
respond to subdividing the 3-cube into three pyramids whose bases are the 
square facets disjoint from one vertex. The six inequalities like Xqoo + 2^001 ^ 4 
correspond to subdividing the 3-cube into two triangular prisms. 

Since the hyperdeterminant -D222 is a factor of the principal determinant 
-E'222, its Newton polytope Af{D222) is a Minkowski summandoi the secondary 
polytope N'{E222)- This implies the existence of a natural map from the 74 
vertices of N{E222) onto the six vertices of A/'(-D222)- Equivalently, we have a 
map from the regular triangulations of the 3-cube to the extreme monomials of 
-D222- A formula for this map can be derived from Theorem 3.2 in §11. A]. 

The extreme monomials of -D222 come in two i?3-orbits. We now describe the 
corresponding two orbits of D -equivalence classes. Each class consists of all 
triangulations which are mapped to a fixed extreme monomial of -D222- The In- 
equivalence class of the monomial cooiCoioCiooCm consists of only one triangu- 
lation, namely the triangulation of type 1 with GKZ vector (1, 5, 5, 1, 5, 1, 1, 5). 
The D-equivalence class of Cqqqc\^^ consists of all 18 triangulations which use 
the diagonal from (0,0,0) to (1,1,1). Thus the number of triangulations of 
the 3-cube decomposes as follows into the sizes of the D-equivalence classes: 

(3) 74 = 2- 1 + 4-18 = 1 + 1 + 18 + 18 + 18 + 18. 

We emphasize that all of these results on the 3-cube are easy and well- 
known. This section served a purely expository purpose, namely, to set the 
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stage for the new results on the 4-cube to be presented in the later sections. 
For instance, Table El generalizes the identity © from the 3-cube to the 4-cube. 

3. SCHLAFLl'S FORMULA AND ITS EXPANSION 

We now consider a multilinear polynomial in four variables x, y, z, and w, 

F = coooo + CoQoiW + Cqoioz + cooiizw + coiool/ + CoiQiyw 
+ coiioyz + coiuyzw + Cioooa; + Cwoixw + Cwioxz 
+ ciouxzw + cnooxy + cnoixyw + cnioxyz + cimxyzw, 

where the 16 coefficients Cijki are regarded as unknowns. 

Definition 3. The hyperdeterminant -D2222 is the unique irreducible polyno- 
mial (up to sign) of content one in the 16 unknowns Qjki that vanishes when- 
ever the system of equations 

(4) F- — - — - — -0 

dx dy dz dw 

has a solution {xo,yo, Zo,Wo) in C^. 

In theory, we can compute -D2222 by eliminating the four variables x, y, z, w 
from the five equations (^J, say by using Grobner bases or resultants J^^, but in 
practice this is infeasible. The analogous computation for a multilinear poly- 
nomial in three variables, however, is easy to do, and it yields the expression 
for D222 stated in the Introduction. 

Schlafii's formula for -D2222 is obtained as follows. We replace each of the 
eight unknowns Ctjk in the 2 x 2 x 2-hyperdeterminant -D222 by Cijko + Cijkiw. 
The resulting expression is a polynomial of degree four in the variable w. Its 
coefficients are polynomials of degree four in the Cijki- The discriminant of this 
polynomial with respect to the variable w is an expression of degree 24 in the 
Cijki- This expression is Schlafii's formula for the hyperdeterminant -D2222- 

Proposition 4. The hyperdeterminant -D2222 coincides with the discriminant 
of D222{cijko + Cijkiw), considered as a polynomial in w, divided by 256. 

This proposition is a special case of jT^i Theorem 14.4.1]. Its proof is based 
on the fact that both polynomials have degree 24, which we know for -D2222 by 
O Corollary 14.2.10]. For tensors of larger format (e.g., 2x2x2x2x2) the 
same recursive method does not work so well. In general, Schlafii's formula 
yields the desired hyperdeterminant times a large extraneous factor. 

The first assertion in Theorem 1 was proved by expanding Schlafii's formula, 
as follows. Using MAPLE we expressed -D222(cjifco + Cijkiw) as 

-D222(cijfco + Cijkiw) = b^w'^ + bsw^ + b2W^ + biw + b^ 
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where each bi is a degree 4 polynomial in the 16 unknowns Cijki- These expres- 
sions for the bi in terms of the Cijki were then substituted into the discriminant 

2566^63 _ I92blbibsbl - 1286^6^^^ ^ lUblbibjb^ - 27blbl 
(5) + lUbobjbibl - Gboblbjb^ - SObobibjbsb^ + 186ofci&2&3 + 1660^2^4 
- Aboblbl - 21b\bl + 186?&2&3&4 - ^h\bl - ^blblb^ + blblbl 

After substituting, MAPLE was unable to expand and combine the discrim- 
inant's 16 products due to memory constraints. Instead, wc used MAPLE to 
substitute and expand each of the 16 products (such as b\b\b1 ) separately. 
We then wrote a C-I--I- program to merge the 16 expansions and divide all 
coefficients by 256. The result of this merge was a lexicographically sorted 
list of 2894276 monomials of degree 24 in the 16 unknowns Cijku each with its 
integer coefficient. The largest absolute value of any coefficient is 112464. 

The symmetry group of the 4-cube is the Weyl group i?4 of order 384. 
It acts on the set of monomials in -D2222 giving 9617 orbits. At our web- 
site, the 9617 lexicographically smallest monomials in each orbit are listed 
in lexicographic order. Each monomial is listed on a separate row, in the 
format [[ExponentVector], Coefficient, FaceDimension, OrbitSize]. Here 
FaceDimension is the dimension of the smallest face of A/'(-D2222) containing 
the exponent vector. For example, the 192 monomials in the orbit of 

2 2 3 5 7 
~2 ■ CooiiCoioiCoiioCoiiiCiQooCloOlClOloClOllCllOoCllOl 

are represented by the row [[O, 0, 0, 2, 0, 2, 3, 5, 7, 1, 1, 1, 1, 1, 0, O], -2, 3, 192]. 
The distribution of the 9617 orbits according to orbit size is as follows: 

2 8 12 16 24 32 48 64 96 192 384 
2 6 6 16 27 24 142 90 577 2743 5984 

The two orbits of size two are represented by the monomial 

2008 • CoooiCooioCoioo'^oiii'^iooo'^ioii'^iioi'^iiiO' 
whose exponent vector lies on a 3-dimensional face, and the monomial 

112464 • cooooCoooiCooio'^ooiiCoiooCoioiCoiioCoiii 
2 222 

■ Ciooo'^lOOlClOloCioiiCllooCiioiCiiio'^llll' 

which has the largest coefficient and whose exponent vector hes in the relative 
interior of A/'(-D2222)- The largest odd coefficient appears in the monomial 

42 2 22222 

— 5811 • CooOlCoOloCooilCoiooCoiOlCoiioCoillCioooCloOlClOloCllOOCllOlCllloCllll 

with orbit size 384 whose exponent vector hes in the relative interior of a 9-face. 
The distribution of the 9617 orbits according to face dimension is as follows: 

0123 4 5 6 7 89 10 11 
111 230 269 540 1145 1862 2138 1845 976 405 70 26 
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The computation of the above face dimensions and the computation in the 
rest of this section were all performed after we computed the facet represen- 
tation of the Newton polytope A/'(-D2222)- In the next section, we present all 
eight classes of facet inequalities and explain how we obtained them. 

In the remainder of this section we discuss the 20992 monomials whose ex- 
ponent vectors lie in the Newton polytope A/'(L'2222) but which have coefficient 
zero in the expansion of hyperdeterminant -D2222- They come in 69 -B4-orbits. 
In the following table we list one representative from each orbit: 



0002114445101010 
0002115345101001 
0002122571111001 
0002216126110110 
0002411450141100 
0003114353100210 
0003144050220201 
0003225033121101 
0004114252100310 
0004215043100202 
0006113141102310 
0011113553111100 
0011222462020011 
0012310542212001 
0013222223212110 
0014320252000032 
0111224132130210 
0121212324202002 



0002114454011010 
0002115354011001 
0002126126200110 
0002223353200002 
0003113444102010 
0003114361110300 
0003222350221200 
0003324041130201 
0004115153000211 
0005114151100410 
0006114060000411 
0011121672100011 
0012125136000021 
0012410441042010 
0013330251110022 
0014330151100023 
0112214132031210 



0002114454100110 
0002115354100101 
0002136025300101 
0002226026200002 
0003114344012010 
0003115262010201 
0003223225201020 
0003401441132100 
0004115161010301 
0005115050110500 
0006212122114110 
0011121681010011 
0012142241311101 
0013143052100013 
0014220351011130 
0014331051010023 
0113213132201202 



0002114463010110 
0002115363010101 
0002144114501100 
0002314235200002 
0003114353012001 
0003134114402100 
0003223242211101 
0004113371000211 
0004125043011111 
0005115060010401 
0011111771111100 
0011144151130101 
0012220551212001 
0013221333201120 
0014231151010131 
0023231152000023 
0114213031112202 



For instance, the underlined vector represents the monomial 

7 7 

CoOloCoOllCoiOOCoiOlCoiloCoiiiCioooClOOlClOloClOllCllooCllOl- 

The orbit of this monomial has size 96. Each of these 69 monomials listed 
above does appear during the expansion of Schlafli's formula, namely, when 
we expand each of the 16 terms of the discriminant (jSJ. But the coefficients 
of these monomials sum to zero when we merge to form -D2222- 

The distribution of the 69 orbits of "missing monomials" according to the 
face dimension in A/'(-D2222) is as follows: 

01 234 56789 10 11 
00 15 3 20 13 7560 

The following method was used to generate all lattice points in A/"(iI'2222) 
and hence find the 20992 "missing monomials" in 69 orbits. We used the five 
equations A ■ x = (24, 12, 12, 12, 12)^ and 268 facet defining inequalities of 

^^{D 

2222) which will be presented in the next section. By symmetry, it suffices 
to generate at least one lattice point in each B4-orbit. We first listed all ways to 
assign eight nonnegative integers summing to 12 on the first facet, i.e., the first 
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eight coordinates. There were 50388 ways in all. Per the symmetry remark 
above, we kept only the 1349 ways which were lex-min under B^. We then 
extended each of these 1349 facet assignments to the entire 4-cube as follows. 
We fixed two disjoint non-parallel edges on the opposite facet, assigned all 
possible values for those four entries, and solved for the remaining four entries 
via the five linear relations above. Of the resulting 4-cube assignments, 87435 
have non-negative entries. We then tested these against the facet inequalities 
and found that 80788 of them lied in the polytope. We now had at least one 
representative from each orbit of lattice points in M {02222)- From these we 
found 9686 orbits of lattice points. Removing the 9617 orbits of terms which 
do appear in 1^2222, we found precisely the 69 additional orbits listed above. 

4. Computing the Newton polytope 

In this section we present our census of the vertices, facets, and other faces 
of the Newton polytope of the 2 x 2 x 2 x 2-hyperdeterminant. The following 
result completes the proof of Theorem 1. Being the analogue of equation ((H), it 
shows the increase in complexity when passing from the 3-cube to the 4-cube. 

Proposition 5. The f -vector of the Newton polytope of D2222 equals 

f{M{D2222)) = (25448, 178780, 555280, 1005946, 1176976, 

927244, 495936, 176604, 39680, 5012, 268). 

We computed the vertices of the Newton polytope A/'(-D2222) starting from 
the exponent vectors of the 2894276 monomials in 9617 orbits described in 
Section 3. We took advantage of the symmetry group using the following 
heuristics: First, we sorted the set of 9617 lex-min elements, one from each 
orbit, lexicographically and then chose a small lexicographically contiguous 
subset. Within this subset, we removed points that were redundant (i.e., not 
vertices) using POLYMAKE ^1]. Whenever a point was found to be redundant, 
we removed its entire i?4-orbit from our original set of points. After a few itera- 
tions we obtained the convex hull of the 9617 lex-min elements, which has 1794 
vertices. The union of the orbits of those 1794 points contains 484804 points. 
We repeated the same process on these 484804 points, finding redundant points 
in a small subset and removing their whole orbits. We eventually reached a 
subset of 25448 points in 111 orbits which appeared to be irredundant. We 
verified that these points are vertices by solving 111 linear programming fea- 
sibility problems. 

The tangent cone at a vertex v is the set {v + w & M}^ : v + e ■ w G 
N'{D2222) for some e > 0}. The faces of the tangent cone are in natural bijec- 
tion with the faces of J\f{D2222) containing v. The vertex figure of a vertex v is 
the 10-dimensional polytope obtained by cutting off v from its tangent cone. 

Proof of Proposition\^ The facets and /-vector of A/'(-D2222) were computed 
as follows: For each of the 111 vertex classes, we chose a representative v and 
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subtracted it from each of the other 25447 vertices. The cone generated by 
these 25447 difference vectors is (a translate of) the tangent cone at v. Using 
POLYMAKE we computed the faces of all dimensions of the 111 tangent cones. 
These computations were fairly fast since A/'(-D2222) is close to being simple. 
We now had at least one representative from the i?4-orbit of every face of 
2222)- Each face is represented by the set of facets containing that face. 
With careful relabeling, we merged the 111 lists of faces and computed their 
orbits, thus obtaining all faces of A/'(D2222)- D 



Here is a complete list of representatives for the 
upper and lower indices attached to each vector 
orbit size respectively. 

0002022662202000.^ 
0002022690010110 
0002023580021100 



111 orbits of vertices. The 
are the coefficient and the 



3^ 

-1 

-'384 



0002022671 1 12000^92 00020226801 1 1 100 

0002023572101001 



0002023571022000^84 



0002024480020200} 
0002027127101010} 



192 



000202358101 lOOlgg^ 0002023590010101 

0002027126112000; 



192 
1 

192 



192 



0002026226202000(^6 
0002027135022000^84 



192 

1 

384 
1 

192 



00020271450101101 
0002028035021100 
0002033470032000 
0002036136100020^ 
0002037026111010 



0002037036100011384 



000202. 



^111100^92 



0002027136011010^84 



192 



384 
1 

384 



0002028036011001^84 

0002033481100002^6 
000203614501002038^4 
-1 

192 



0002037034032000 



0002037045010011^84 



192 



0002046060040200 }f:o 
0002111771112000^6 



0002111771201100}92 



0002028027101001 _ 

0002028045010101^92 
000203349001000296^ 
0002037025122000^92 
0002037035021010^84 
0002037045100002 
0002037054010002^9^2 0002046035110020384 0002046044020020384 

0002055050052000 o«^^ 0002055070030002 o6^^ 

0002111780111100384 
0002112671022000384 

0002112680021100^84 
0002112681100101^84 
0002117117201100}92 
0002117126200200384 
0002117136010110384 
0002118027100101384 
000212258110000296^ 

'192 
'192 



000211178110011070^0 0002111790010110} 



•"192 

000211267120020038^4 

0002112680110200^84 
-1 

192 
1 

384 



192 

-1 
384 



192 



0002112690010101 
0002117118101010 



0002117127011010^84 
0002118018101001}qo 0002118027011001 



0002112672101001 
0002112681011001} 
0002117117112000^ 
0002117126022000 



0002117127100110^84 



96 
1 

192 



'-192 

0002118036010101^84 
0002122590010002^8 
0002201772100110^84 



384 

0002122570032000}92 
0002201763101010384 
0002202662022000^6 



0002201772011010 
0002202662200200 
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0002202663101001^84 



0002217036010002 
0002300762200110 
0002302563100002 



-1 

192 
-1 
192 
-1 
192 



000220267210010138^ 

0002300752301100^84 
0002300771110110^92 
0002306136010002^6 



0002217027100002^ 



96 



0002300761211100^84 
000230166310001li 



384 



0002400652300101^84 



000240066121010l}92 0002433050050002q«i 0002500552210002,^ 



96 



-48 



0003033390000003^ 
0003036045000021 
0003111671102001 

0003117017101110^ 



16 
1 

192 
-1 
192 



192 
1 
64 



96 



0004044040044000^ 
0009111011107110^ 
001111177210100lj 



96 



48 



00030351450000301 
0003036063000003 
00031116801011011 
0003117026100210 



0004044080000004 ]^^ 



96 



192 
1 

192 



0003036026101020^ 
0003045035100030}92 
000311169000011164^ 

000311703600011119^2 



0008111111118000^ 
0011110881011010^ 
0011120771201001} 

0011120781100011^9^2 0110100990010110^ 0110106336010110^ 



48 



192 



0009111011108001 
0011111781011001 



1 

96 
-1 
192 



92 



192 



32 



192 

44 



384 
39 



The distribution of the 111 orbit sizes is 

8 16 32 48 64 96 
1 2 2 4 3 16 

The distribution of the 111 corresponding extreme monomial coefficients is 

16 -16 -27 
1 2 1 



1 

60 



-1 
47 



All but four of the monomials have coefficient ±1. The largest absolute value 
of a coefficient appears in the monomial — 27-CQoiiCQiQiCgiiQcfQQQCim which has 
orbit size 16. Its D-equivalence class will be discussed in Example E| The 
monomial 16 ■ CqohCoioiCohoCioooCiohCiioi has orbit size 192, and there are two 
monomials with coefficient —16, each having orbit size 96. Their D-equivalence 
classes will be discussed in Example UHl 

Among the 111 types of vertices, 35 are simple. The last underlined ver- 
tex, representing the monomial 8 ■ CoooiCooioCoiooCoinC^oooCioiiCiioiCino, has the 
largest vertex figure. That vertex figure has the /-vector 

(67, 873, 4405, 11451, 17440, 16452, 9699, 3446, 667, 56). 

Thus this distinguished vertex is adjacent to 67 other vertices, and it lies in 56 
facets. In Example ITHl we shall see that this vertex corresponds to the largest 
D-equivalence class of triangulations. Table 1 shows the distribution of the 
111 types of vertices according to the number of incident edges and facets. 

The 268 facets of M{D2222) come in eight orbits. In Table 2 we present 
the /-vector and the orbit size for each facet. We discuss the geometry of 
each of the eight facet types. Since M{D2222) is a Minkowski summand of 
the secondary polytope J\f {E2222), each facet normal of A/'(-D2222) is a facet 
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total 


35 


28 


11 
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1 
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111 



Table 1. Distribution of the 111 vertices according to the num- 
ber of facets (columns) and edges (rows) containing the vertex. 





/-vector 


orbit 


1 


(11625,72614,197704,308238,303068,194347,80874,20906,3021,187) 


16 


2 


(4112, 25746, 71456,115356,119228,81590,36802,10488, 1704, 122) 


12 


3 


(363. 2289. 6538. 10996. 11921. 8581. 1080. 1239. 225. 22) 


8 


4 


(938, 5226, 13182, 19716, 19263, 12765, 5758, 1721, 318, 31) 


32 


5 


(336, 1937, 5126, 8121, 8468, 6022, 2928, 950, 194, 22) 


48 


6 


(289, 1624, 4228, 6636, 6894, 4914, 2413, 798, 168, 20) 


96 


7 


(450. 2526. 6522. 10103. 10315. 7195. 3440. 1099. 220. 24) 


48 


8 


(681. 3906. 10323. 16107. 17191. 12261. 5933. 1877. 357. 31) 


8 



Table 2. The eight types of facets of the Newton polytope A/'(-D2222) 



normal of A/'(-E2222) and hence corresponds to a coarsest regular subdivision of 
the 4-cube. We depict each subdivision by its tight span, which is the complex 
of bounded faces of the polyhedron {u E M.^ : u ■ A < w}, where A is the 
5 X 16-matrix in Section 1 and w e M}^ is the facet normal in question. 

Facet 1: These 16 facets are defined by inequahties hke xoooo > 0. The tight 
span of the corresponding subdivision is a line segment. The two maximal 
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cells are the simplex obtained by slicing off one vertex and the convex hull of 
the other 15 vertices, which has /-vector (15, 34, 28, 9). 

Facet 2: These 12 facets are given by inequalities hke 

2^0000 + 3:^0001 + 2^0010 + 2:^0011 ^ 2. 

The tight span is a line segment. The two maximal cells in the corresponding 
subdivision have /-vector (12,24,19,7). Here we cut the 4-cube with the 
hyperplane containing the eight vertices in a pair of opposite 2-faces. 

Facet 3: These 8 facets are given by inequahties like 

2 • 3^0000 + 3:^0001 + 2^0010 + 3:^0100 + 3:^1000 ^ 2. 

The tight span is a line segment, and the /-vector of both maximal cells is 
(11, 28, 26, 9). Here we cut the 4-cube with the hyperplane spanned by the six 
vertices which have coordinate sum two, whose convex hull is an octahedron. 

Facet 4: These 32 facets are given by inequalities like Xqooo + s^oooi ^ 9. The 
tight span is a triangle. Each of the three maximal cells is a prism over a 
square pyramid and has /-vector (10, 21, 18, 7). They are formed by fixing an 
edge of the 4-cube and joining that edge to the three facets which are disjoint 
from that edge. Any two cells intersect in a triangular prism, and all three 
cells intersect in the rectangle formed by the fixed edge and its opposite edge. 

Facet 5: These 48 facets are given by inequalities hke 

3^0000 + 2^0001 + 2^0010 + 3:^0100 + 3:^0110 + 3:^1000 + 2^1001 ^ 3. 

The tight span is a triangle. Two of the three maximal cells have /-vector 
(10, 23, 21, 8) and intersect in a square pyramid; the other cell has /-vector 
(9, 18, 15, 6) and intersects each of the first two in a triangular prism. The 
intersection of all three cells is a square. 

Facet 6: These 96 facets are given by inequalities like 

3^0000 + 2^0001 + 2:0010 + 2:0011 + 2:0100 + 2:0110 + 2:1000 + 2:1001 ^ 4. 

The tight span is a tetrahedron. Two of the maximal cells have /-vector 
(8, 18, 17, 7) and the other two have /-vector (9, 21, 20, 8). The intersection of 
all four cells is a segment which is a diagonal of the 4-cube. 

Facet 7: These 48 facets are given by inequalities like 

2:0000 + 2:0001 + 2:0010 + 2:0011 + 2:0100 + 2:1000 + 2:1100 ^ 19. 

The tight span is a square pyramid. The maximal cell corresponding to the tip 
of the pyramid is the convex hull of two squares in complementary dimensions 
that share a vertex. It has /-vector (7, 17, 18, 8) and meets each of the other 
four cells in a triangle. The other four maximal cells have /-vector (8, 18, 17, 7), 
and their intersection is a triangle. All the cells meet along a diagonal segment. 
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(1,5) (23,15) (12,12) (12,12) (12,11) (12,11) 
• • • • • • 

facet 1 facet 2 facet 3 




Figure 2. The tight spans dual to the subdivisions of the 4- 
cube corresponding to the facets oi N {02222)- Each vertex of a 
tight span is labeled by the normalized volume and the number 
of vertices of the corresponding maximal cell in the subdivision. 



Facet 8: These 8 facets are given by inequahties like 

2 • 3^0000 + 2^0001 + 3^0010 + 3^0100 + 2^1000 ^ 18- 

These facets are indexed by diagonals. We can identify the 4-cube with a 
boolean lattice when wc fix a diagonal. The maximal antichain in this lattice 
has the structure of an octahedron. Each point in the antichain uniquely de- 
termines a maximal cell in the subdivision. Hence the tight span is also an 
octahedron. Each maximal cell is the convex hull of two squares in comple- 
mentary dimensions joined at a vertex, and each has /-vector (7, 17, 18, 8). 

This concludes our discussion of the facets of the Newton polytope M {02222)- 
The i?4-orbits of the eight stated inequalities yield 268 facet inequalities. These 
together with the identity A ■ x = (24, 12, 12, 12, 12)^ gives an irredundant 
presentation of N{D2222) by linear equations and inequalities. 
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5. Tight spans and A-discriminants of other 0/1-polytopes 

The concepts and computations presented in this paper make sense for any 
d X n-integer matrix A whose row span contains the all-ones vector (1,1,. ..,1). 
Following ^21; the A-discriminant Da is an irreducible factor of the principal 
A-determinant Ea- The vertices of the secondary polytope M^Ea) correspond 
to regular triangulations of A, and these map onto the vertices of the Newton 
polytope M^Da)- The fibers of that map (i.e., the D-equivalence classes) and 
the normal fan of N'{Da) have recently been characterized in [S]. 

For fixed A and an arbitrary row vector w G M", we consider the polyhedron 

= {u eR'^ : u- A<w}. 

This polyhedron is dual to the regular triangulation of A defined by w. The 
tight span of {A, w) is the complex of bounded faces of the polyhedron P^. For 
experts in tropical geometry, we note that the [d — 2)-skeleton of P^ is the 
tropical hypersurface of {A, w). The polyhedron P^ also occurs naturally when 
relating regular triangulations to linear programming duality as in pUl §1.2]. 

We propose that, for many matrices A of interest in combinatorics and 
discrete convexity the tight span is an excellent geometric representation 
of the data {A, w). In Section 2 we studied the case when A is the 4 x 8-matrix 
representing the 3-cube, and in Sections 4, 6 and 7 we are concerned with the 
case when A is the 5 x 16-matrix representing the 4-cube. In what follows we 
present two related situations which have appeared in the recent literature. 

Example 6 (Tropical Polytopes). Let A be the (r + s) x (r ■s)-matrix which 
represents the direct product of an (r — l)-simplex A^-i with an (s — l)-simplex 
As_i. The principal A-determinant Ea is the product of all subdeterminants 
(of all sizes) of an r x s-matrix of unknowns (cjj); see ^3 page 303]. 

For instance, if r = s = 3 then A represents the product of two triangles, 

■ 



1 1 

' 

1 
1 . 

and the A-discriminant Da is the determinant of the matrix 

Cll Ci2 Ci3 
C21 C22 C23 
C3I C32 C33 

The degree 30 polynomial Ea is the product of all 19 subdeterminants of this 
matrix. The extreme monomials of Ea appear in ^21 Figure 39, page 250]. 

The Newton polytope N{Ea) of the product of all subdeterminants is the 
secondary polytope whose faces correspond to regular polyhedral subdivisions 
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of A.f._i X As_i. The corresponding tight spans are precisely the tropical poly- 
topes, which are obtained as tropical convex hulls P of r points in (s — l)-space. 
The combinatorial classification of such polytopes is the analogue of what will 
be accomplished for the 4 cube in the next section. We refer to recent work of 
Santos IIH] for a discussion. His Section 5 deals with tropical polytopes, and 
his Figure 2 is essentially the same as our Figure 1. Our number 235277 in 
Theorem 2 is the 4-cube analogue of the numbers in the table at the end of 
Section 4 in including the number 35 in Fig. 1] and in ^ Fig. 6]. □ 

Example 7 (Injective hulls of finite metric spaces). The term "tight 
span" is derived from the special case when A represents the second hypersim- 
plex. As shown in [1^120]! here the tight spans are precisely the injective hulls 
of finite metric spaces, which play an important role in phylogenetic combina- 
torics. The case studied in [201 concerns metrics on six points, where 

■111110000000000- 
100001111000000 
. _ 010001000111000 

001000100100110 ■ 

000100010010101 
.000010001001011. 

Here the A-discriminant is the determinant of the symmetric matrix 

/ Ci2 Ci3 Ci4 Ci5 CigX 

Cl2 C23 C24 C25 C26 

Cl3 C23 C34 C35 C36 

Ci4 C24 C34 C45 C46 

Cl5 C25 C35 C45 C56 

\Cl6 C26 C36 C46 C56 / 

Up to a constant, the principal A-determinant Ea_ is the product of all principal 
minors of size > 3. This is a polynomial of degree 156 having 194160 extreme 
terms in 339 symmetry classes ,20^ Theorem 1]. Thus our number 235277 is 
the 4-cube analogue to the number 339 of generic six-point metrics. □ 

The prominent role of the second hypersimplex in phylogenetic combina- 
torics raises the question of what happens for other hypersimplices. We propose 
the following problem for further mathematical and computational research. 

Problem 8 (The hypersimplex A(6,3)). 

Let A be the 6 x 20-matrix whose columns are the vertices Ci + ej + Ck (1 < 
i < j < k < Q) of the third hypersimplex A(6,3). What is the degree of the 
y4- discriminant Da ? Can the 14-dimensional polytope M{Da) be computed ? 
Can the monomial expansion of the A- discriminant be computed ? 

All matrices A in the examples above have their entries in {0, 1}, so they 
represent subpolytopes of cubes of appropriate dimensions. In what follows we 
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examine the A-discriminants of various subpolytopes of tlie 4-cube. Tliese arise 
naturally as irreducible factors in the initial forms of the hyperdeterminant. 

For any vector w in M^^ we define the initial form in^„(L)2222) as the sum of 
all terms in -D2222 having minimal w-weight. If w is generic then in^(D2222) 
is a monomial, and we have classified all of these monomials in Section 4. We 
now consider the other extreme case when w is as non-generic as possible. 
More precisely, we pick w among the normal vectors to the facets of A/'(-D2222) ■ 

Proposition 9 (All mciximal initial forms of the hyperdeterminant). 

The following list is the classification of the initial forms of hyperdeterminant 
-D2222 corresponding to all facets of its Newton polytope. The eight symmetry 
classes of facets of J\f {D2222) are listed in the same order as in Section 4- 

Facet 1: If w = eoooo then in^ (1^2222) is obtained from £^2222 by setting 
coooo = 0. This initial form is irreducible: It is the 74-discriminant where A is 
the 15-point configuration obtained from the 4-cube by removing one vertex. 

Facet 2: If w = eoooo +eoooi+eooio + eooii then in^(£)2222) has 67230 terms and 

factors as CioiqCiooi — ciooqCioh times CqioiCoho — CoiooCoiii times the square of 
the 2 X 2 X 2-hypcrdctcrminant of [ciqoo, Ciqoi, Cioio, Cion, cqiqo, cqioi, Cquo, cqih] 
times the product of two larger factors, each having 66 terms of degree six, 
which are the 74-discriminants of the two maximal cells in the subdivision. 

Facet 3: If w = 2 • eoooo + eoooi + eooio + eoioo + eiooo then in^(i:)2222) equals 



(1/4) • cooiicoioicoiiociooicioiociioo • { det 

2ciiii ciiio ciioi cioii Com 

ciiio ciioo cioio Cquo 

CllOl Ciioo Ciooi CoiOl 

ClOll ClOlO ClOOl Cooil 

Com coiio coioi cooii 



Cooil CoiOl CloOl 

COOII CoilO ClOlo 

coioi coiio ciioo 

ClOOl ClOlO Ciioo 



det 



2coooo Coooi CooiO Coloo Clooo 

CoOOl Cooil Coioi Clool 

coolo Cooil coiio Clolo det 

COlOO coioi coiio ciioo 

ClOOO CloOl ClOlO Ciioo 



The squared factor is the 74-discriminant of the octahedron, and the 5x5- 
determinants are A-discriminants of the two maximal cells in this subdivision. 
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Facet 4: li w = -eoooo - eoooi then in^(-D2222) equals 

(coooociiii — CooOlClllo)^ 

(cilOOCllll — CiioiCiiio)(cioioCllll — CioilCiiio)(coil0Cllll — CoillCiiio) 

('^oooo'^iooiciiii — CqoooCioiiChoi — coooocoooiciooociiii — coooocoooiciooiciiio 

+COOOOCOOOICIOIOCIIOI + CooOOCOOOlClOllCllOO + CoooiClOOOClllO — CoooiClOlOCllOOj 

('^OOOO'^OlOlCllll — CooooCOlllCllOl — coooocoooicoioociiii — coooocoooicoioiciiio 

2 2 \ 

+COOOOCOOOICOIIOCIIOI + CooOOCOOOlColllCllOO + CqooiCOIOOCIHO — Cqqqj^CoiioCiiooJ 

/ 2 2 

I'^oooo '^0011 ci 111 — CqqqqCoiiicioii — coooocoooicooiociiii — coooocoooicooiiciiio 

2 2 \ 

+C0000C0001C0110C1011 + CooOOCOOOlCoillClOlO + CqooiCooiociho — CoooiCoiiocioioj- 
The quartic factors are the 74-discriminants of the three maximal cells. 

Facet 5: If w = eoooo + eoooi + eooio + eoioo + eoiio + eiooo + eiooi then in^(L'2222) 
equals cooooCooiiCoioiCioioCnoo times (coonCnoo - coioiCioio)^ times 

CoOllCllOOCllll — CooilCiioiCiiio — CoioiCioioCllll 
+C0IOIC10IIC1110 + CoillCioioClloi — CoillCioilCiloo 

times the product of two larger factors, each having 15 terms of degree five 
which are the 74-discriminants of the two cells with /-vector (10, 23, 21, 8). 

In the remaining cases the tight spans are three-dimensional (Figure 2), and 
the initial form in^(D2222) factors as a monomial times the product of the 
74-discriminants of the maximal cells in that subdivision of the 4-cube. 

Facet 6: If w = eoooo + eoooi + eooio + eooii + eoioo + eoiio + eiooo + eiooi then 
in^ (1^2222) has 404 terms and factors as a monomial times the product of the 
^-discriminants of the four maximal cells: 

in^(D2222) = Coooo ■ Ciioo • cfoio • Cqioi ■ 

(coiOlClOloCllll — CoioiCioilCiiio — CoillCioioClloi + CoillCioilCiloo) 

(coOOOCoillCioil — CoooiCoillCiolO ~ CooioCoiOlClOll + CooilCoiOlClOlo) 

(coOOOCoiOlCoillCiiio — CooooCoiiiCiloO ~ CooioCoioiClllO 

2 \ 

+C0OIOC0IOIC0IIIC1100 ~ CoiooCoiOlCoillCiolO + CoioiCoiloClOloj 
(coOOOClOloClOllClloi — CooooCioiiCllOO ~ CoooiCioio^llOl 

2 \ 

+C0OOIC10IOC10IIC110O ~ CoioiCioooClOloClOll + CoiOlCloOlCioioj- 
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Facet 7: li w = -eoooo - eoooi - eooio - eooii - eoioo - eiooo - enoo then 

in^(D2222) = ciiii ■ Cqoh ■ cfioo ■ 

(cooooCooiiCiioo ~ CoooiCooioCllOO ~ CooilCoiooClOOo) 
(cooioClOOOCllll — CooioClOllCllOO ~ CooilCioooClllO + CooilCioioCllOo) 
(cooioCoiOOCllll — CooioCoillCiloO ~ CooilCoiooClllO + CooilCoiloCllOo) 
(CooOlClOOOCllll — CoooiCloilCiloO ~ CooilCioooCllOl + CooilCiooiCllOo) 
(cooOlCoiOOCllll — CooOlCoillCiloO ~ CooilCoiooCllOl + CooilCoiOlCllOo)- 

Facet 8: liw = -2- eoooo - eoooi - eooio - eoioo - eiooo then in^(D2222) equals 

CqoOO ■ '^1111 ■ (cooooCiiooCllll — CooOOCiioiCiiio — CoiooClOOOCllll) 
(CoOOOClOloCllll — CooooCloilCiiio — CooioCloOoCiiii) 
(coOOOCoiloCiiii — CooooCoillCiiio — CooioCoiooCiiii) 
(CoOOoCiooiCiiii — CooooCloilCiioi — CoooiCloOoCllll) 
(CoOOoCoiOlCllll — CooooCoillCiioi — CoooiCoiOOCllll) 
(coOOOCooilCiiii — CooOOCoillCioil — CoooiCooioCiiii). 

These formulas show that A- discriminants for sub configurations of the cube 
appear naturally as irreducible factors of leading forms of the hyperdetermi- 
nant. This suggests the general problem of studying A-discriminants and the 
related tight spans for various families of 0/1-polytopes. Other natural classes 
of configurations for which a study of A-discriminants would be interesting in- 
clude generalized permutohedra, Birkhoff polytopes, and reflexive polytopes. 



6. Vertices of the secondary polytope 

The 4-cube is given by the 5 x 16-matrix A in the Introduction. For any 
generic vector w G M}^ we consider the tight span of {A, w) whose cells are the 
bounded faces of the simple polyhedron = {u ^ MP : u ■ A < w}. Each 
vertex of is indexed by five columns of A. This collection of 5-tuples, each 
regarded as a 4-simplex, is a regular triangulation 11^ of the 4-cube (cf. fOl 
§5]). 

The regular triangulations of the 4-cube are in bijection with the vertices 
of the secondary polytope. As in Section 2, we define the secondary polytope 
as the Newton polytope A/'(-E'2222) of the principal determinant E22221 which is 
the following product of determinants associated to all faces of the 4-cube: 

• 16 linear factors Cijki corresponding to the vertices, 

• 24 quadratic factors D22 corresponding to the 2-faces, 

• 8 quartic factors 1^222 corresponding to the facets, 

• the hyperdeterminant -D2222 corresponding to the solid 4-cube. 
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Figure 3. Tight spans of the three triangulations in Example ITTl 

Thus the principal determinant -E2222 is a polynomial in the 16 unknowns Cijki- 
of degree 120 = 16 ■ 1 + 24 ■ 2 + 8 ■ 4 + 1 ■ 24. The vertices of its Newton polytope 
7V(-E'2222) are the GKZ vectors of the regular triangulations of the 4-cube. 

Recall from the first part of Theorem 2 that the 4-cube has 87959448 regular 
triangulations, or, equivalently, ^{£2222) has 87959448 vertices. Up to the 
action of the Weyl group -B4, there are 235277 types of regular triangulations. 
Thus, this is the number of combinatorially distinct simple polyhedra Pw, and 
also the number of generic tight spans. 

The main objective of this section is to provide a detailed study of the 
regular triangulations of the 4-cube, with emphasis on their tight spans. We 
present a range of results which were computed from our data. 

Corollary 10. The distribution of the 235277 symmetry classes of vertices of 
the secondary polytope M{E2222) according to orbit size is as follows: 

8 16 24 32 48 64 96 192 384 
3 7 2 13 48 102 516 11357 223229 

Example 11 (The most symmetric triangulations). One of the triangu- 
lations with orbit size 8 is the staircase triangulation. Its tight span is a solid 
permutohedron, and its GKZ vector is (4, 6, 6, 24, 6, 4, 4, 6, 6, 4, 4, 6, 24, 6, 6, 4). 
Another is the unique triangulation that uses 16 maximal simplices. Eight 
simplices have volume one and eight have volume two. Its GKZ vector is 
(1, 12, 12, 1, 12, 1, 1, 20, 20, 1, 1, 12, 1, 12, 12, 1), and its tight span is a 3-cube 
with tentacles attached to its vertices. The third triangulation with orbit 
size 8 has a sohd truncated 3-cube as its tight span. Its GKZ vector is 
(3,8,8,3,8,3,3,24,24,3,3,8,3,8,8,3). The tight spans of these symmetric 
triangulations are shown in Figure |21 □ 

Table 3 classifies the regular triangulations according to the number of max- 
imal simplices used. It is well-known that the minimal number is sixteen 
(cf. |3]). We see that any triangulation can only use at most one simplex 
of volume three, since all such simplices contain the centroid of the cube as 
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17 
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344 


21 
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21 


19 
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1263 


22 


20 


2 





24499 


22 


21 





1 


1967 


23 


22 


1 





48648 


24 


24 








151989 



Table 3. Classification of regular triangulations of the 4-cube 
according to the number and volumes of maximal simplices 



an interior point. There are 16 simplices of volume three, all of which are 
54-equivalent; a representative is {(1000), (1111), (0011), (0101), (0110)}. 

Example 12 (Triangulating the 4-cube with 17 simplices). Table 3 
shows that there is a unique regular triangulation with 17 maximal simplices, 
ten of volume one and seven of volume two. Its tight span, shown in Figure |3J 
is a cube with a truncated vertex and tentacles attached to the seven original 
vertices. Its GKZ vector is (1, 11, 12, 1, 12, 1, 1, 21, 20, 3, 1, 11, 1, 11, 12, 1). □ 

Table E] classifies the tight spans according to their /-vectors and their sig- 
nature, by which we mean the set of dimensions of the maximal cells. Note 
that every triangulation of the 4-cube has at most one interior edge (namely, a 
diagonal), so each tight span has at most one 3-dimensional cell. The /-vector 
of a tight span is determined by the number of vertices and 3-cells, i.e., the 
number of maximal simplices and diagonals used in the triangulation. 

Example 13 (The smallest tight spans of signature {3, 2} and {3, 2, 1}). 

The unique tight span of signature {3, 2} with 20 vertices contains a solid tetra- 
hedron with six hexagons attached, one to each of the tetrahedron's edges. See 
FigureH Its triangulation has GKZ vector (5, 2, 2, 15, 2, 15, 15, 4, 24, 5, 5, 2, 5, 2, 
2, 15) and lies in the 88th D-equivalence class. 

There are two tight spans of signature {3,2,1} with 18 vertices. Their 
GKZ vectors are (1, 9, 13, 2, 12, 1, 1, 21, 19, 6, 1, 9, 1, 11, 12, 1) and (1, 12, 12, 1, 
11, 1, 1, 21, 19, 1, 1, 13, 5, 10, 10, 1). They are both in the largest D-equivalence 
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{3} 
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(17, 22, 7, 1) 






1 








1 


(18, 24, 8, 1) 
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(19, 26, 9, 1) 
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19 




48 
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(20, 28, 10, 1) 
(20, 27, 8, 0) 
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739 


(21, 30, 11, 1) 
(21, 29, 9, 0) 
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392 
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1411 
5699 


(22, 32, 12, 1) 
(22, 31, 10, 0) 
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84 


9342 


2156 


14772 


9538 
16928 


(23, 34, 13, 1) 




2116 


100 


46432 






48648 


(24, 36, 14, 1) 


125 


27054 




124810 






151989 


total 


125 


29288 


253 


182186 


2572 


20853 


235277 



Table 4. Distribution of the 235277 types of triangulations. 
The rows correspond to the /-vectors and the columns corre- 
spond to the dimensions of maximal faces of the tight span. 



Figure 4. The triangulations in Examples IT^ and IT^ 

class. One of them is shown in Figure |3] The other has a triangle and a 
tentacle in place of the rectangular maximal 2-face. □ 

The first column in TablelUis especially interesting as it ties in with an active 
area of research in combinatorial commutative algebra. Using the language 
introduced in the c?-cube is a Gorenstein polytope, and the tight spans of 
signature {d — 1} are precisely its Gorenstein triangulations. The 3-cube has 
only one Gorenstein triangulation, corresponding to the last type in Figure ^ 

Example 14 (The Gorenstein triangulations of the 4-cube). The 4- 

cube has precisely 125 types of Gorenstein triangulations. Two of them are 
pictured in Figure El The Stanley- Reisner ideal of a Gorenstein triangulation 
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Table 5. Decomposition 
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into £)-equivalence classes 



corresponds to an initial ideal of the Segre variety P-*^ x P-"^ x P-"^ x c P^^ that 
is squarefree and Gorenstein. The minimal free resolution of the Alexander 
dual of this monomial ideal is cellular and is supported on the tight span itself. 
For tropical polytopes (Example IHI) this was shown in i2|, but the same holds 
for unimodular regular triangulations of arbitrary polytopes, including the 4- 
cube. All 125 Gorenstein tight spans are deformed permutohedra, and it would 
be interesting to study their combinatorics in arbitrary dimensions d > 4. □ 

We now return to our discussion of the hyperdeterminant -D2222- Since -D2222 
is a factor of the principal determinant E22221 there is a natural many-to-one 
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map from the vertices of A/'(-E'2222) to the vertices of M {02222)- Two regu- 
lar triangulations are D-equivalent if they are mapped to the same vertex of 
2222)- For the 3-cube the D-equivalence classes were described in Sec- 
tion 2. Table El shows the decomposition of the number 87959448 of regular 
triangulations of the 4-cube into the cardinalities of the orbits of D-equivalence 
classes. Thus the identity expressed in TableElis the 4-cube version of the iden- 
tity for the 3-cube. The 111 orbits of D-equivalence classes are hsted in 
the order of the corresponding orbits of vertices of A/'(-D2222) in Section 4. 

In the triangulations of the 4-cube, the centroid (1, |, |, |, |) is in the rel- 
ative interior of one of the following three kinds of simplices: a 4-simplex of 
normalized volume 3, a tetrahedron of normalized volume 2, or a diagonal. 

Example 15 (The "coefRcient —27 class"). The triangulations which use 
a fixed simplex of normalized volume three form a single D-equivalence class of 
cardinality 82832. This class corresponds to the vertex 0004044080000004^6^^ 
of Af {02222)- The corresponding monomial in 1^2222 has the largest absolute 
value of any coefficient, namely —27, among all the extreme monomials. □ 

Example 16 (The "coefRcient ±16 classes"). The triangulations in which 
the centroid is in a fixed tetrahedron of volume two form a D-equivalence 
class corresponding to one of the extreme monomials of -D2222 with coefficient 
±16, which come in three symmetry classes underlined in Table and on 
page ^2 There are 24 such tetrahedra, all lying in a single i?4 orbit. The 
hyperplane spanned by the special tetrahedron contains eight vertices of the 4- 
cube, and there are four vertices on each side. In each triangulation, one vertex 
from each side is joined to the special tetrahedron. The Hamming distance 
between these two vertices can be two, three, or four, and they correspond 
precisely to the three D-equi valence classes with coefficient ±16. The first 
underlined class corresponds to Hamming distance three where there are two 
choices on one side after the vertex on the other side is chosen. Hence it 
has orbit size 24 x 4 x 2 = 192. The second and third underlined classes 
correspond to Hamming distances four and two respectively and have orbit 
size 24 X 4 = 96. □ 

The other D-equivalence classes contain only triangulations using a diagonal 
of the 4-cube, and all those classes contain some unimodular triangulations. 

Example 17 (The smallest D-equivalence class). This class has only 64 
triangulations and corresponds to the last summand in Table El The tight 
span of each of these 64 triangulations contains a hexagonal prism and six 
lower-dimensional pieces attached to alternating edges of the prism, as shown 
on the left in Figure El Each of these six lower-dimensional pieces is either a 
rectangle or a triangle with a tentacle, whence the number 2^ = 64. □ 

Example 18 (The largest D-equivalence class). This class has 349555 
triangulations and corresponds to the vertex of A/'(-D2222) which has the largest 
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Figure 5. Tight spans of triangulations in Examples ITTf fTHl and fTTH 



vertex figure and the smallest orbit size. It contains all the triangulations in 
Examples imini and m □ 

Using our data, it is easy to study other questions concerning regular trian- 
gulations of the 4-cube. Answering a question which was left open in ^^1; we 
found that the 4-cube does not possess any delightful triangulations. Such a 
triangulation would correspond to a square-free initial ideal of the first secant 
variety of P-*^ x P-*^ x x c P"*^^, which has dimension 9 and degree 64. For 
every regular triangulation 11 we computed the number of 10-element subsets 
of the vertices of the 4-cube which are the union of vertices of two simplices in 
n. That number is bounded above by 64, and 11 would be delightful if equality 
holds. 

Example 19 (The most delightful triangulation). We found that the 
maximum number of 10-sets of vertices which are the unions of vertices of 
two simplices in a regular triangulation is 56. This number is attained by two 
triangulations which are both in D-equivalence classes of type 11. Their GKZ 
vectors are (1, 4, 4, 11, 6, 13, 9, 12, 12, 9, 13, 6, 11, 4, 4, 1) and (1, 4, 4, 11, 6, 13, 9, 
12, 12, 11, 13, 4, 9, 4, 6, 1). Both of the corresponding tight spans have /-vector 
(24, 36, 14, 1) and consist of two maximal edges, eight maximal 2-faces, and 
one 3-cube. The tight span of the first triangulation is shown in Figure |H1 □ 

7. Facets of the secondary polytope 

The facets of the secondary polytope correspond to the coarsest regular 
polyhedral subdivisions 11. Here coarsest means that 11 refines no other proper 
subdivisions. We computed all such subdivisions for the 4-cube. This re- 
sult gives an irredundant inequality representation of the secondary polytope 
N'{E2222) of the 4-cube, analogous to that of the 3-cube in expression (j2I). 

Since the Newton polytope A/'(-D2222) is a summand of the secondary poly- 
tope A/'(-E'2222), all facet normals of A/'(-D2222) are also facet normals oi N {E2222) ■ 
Thus we already know 268 facets of the secondary polytope from Section 4. 
Furthermore the linear equations Ax = (120, 60, 60, 60, 60)"^ hold on Af{E2222)- 
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In order to compute the facets of the secondary polytope of the 4-cube, 
we first wrote a program to compute the GKZ vector of each regular trian- 
gulation of the 4-cube. The information about the triangulations themselves 
was obtained using TOPCOM. We then used the following method to com- 
pute the facets of the secondary polytope, which can be applied to any point 
configuration with symmetry. 

Since J\f {E2222) is an 11-dimensional polytope with 87959448 vertices, com- 
puting its facets is currently out of scope of any general polyhedral software. 
Thus, to compute the facets, we took advantage of the relationship between 
the secondary polytope and subdivisions of the 4-cube. Namely, we know that 
each vertex of A/'(-E'2222) corresponds to a regular triangulation of the 4-cube, 
and, by [121 §7.2.C], each edge {v,w) of N{E2222) corresponds to the circuit 
connecting the two triangulations represented by the vertices v and w. 

Note that, given any triangulation of the 4-cube, there are relatively few 
circuits which contain a face of the triangulation. As a result, each tangent 
cone of A/'(£^2222) is generated by relatively few vectors. Hence we computed 
the facets of A/'(-E2222) by equivalently computing the facets of its tangent 
cones. Although there are 87959448 tangent cones, they come in 235227 orbits 
under the -B4-action. Thus, we only computed the facets of 235227 tangent 
cones; all other facets of J\f {E2222) were obtained by applying the i?4- action. 

We found that Af{E2222) has 80876 facets in 334 orbits. Thus the 4-cube 
admits exactly 80876 coarsest regular subdivisions. The distribution of the 
types of coarsest regular subdivisions according to orbit size is 

8 12 16 24 32 48 64 96 192 384 
2 1 4 2 4 14 16 26 132 133 

Example 20 (The most symmetric coarsest regular subdivisions). 

The three most symmetric coarsest regular subdivisions come from facets of 

^^{D 

2222)- The two coarsest subdivisions with orbit size 8 correspond to Facet 
3 and Facet 8 in Sections 4 and 5. Their tight spans are a segment and an 
octahedron respectively. The unique subdivision with orbit size 12 corresponds 
to Facet 2 in Sections 4 and 5. Its tight span is also a line segment. □ 

The distribution of the orbits according to the number of maximal cells in the 
subdivision and the dimensions of maximal faces of the tight span are shown 
in Table IHl A coarsest subdivision of the 4-cube can have up to 13 maximal 
cells. Here it is no longer the case that the number of vertices and 3-faces 
determine the /-vector of the tight span. There are 33 types of /-vectors. 

Example 21 (The missing split). The 4-cube admits 4 splits, i.e., subdivi- 
sions obtained by slicing the 4-cube with a hyperplane, whose tight spans are 
hne segments. Three of them correspond to Facets 1, 2, and 3 of M {02222)- 
The remaining split subdivides the 4-cube into 2 cells with 8 and 14 vertices 
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Table 6. Distribution of the 334 types of coarsest regular sub- 
divisions. The rows correspond to the number of maximal cells 
and the columns correspond to the signatures of the tight span. 



(4,8) (20,14) 
• • 



Figure 6. The tight spans for Examples EH E21 and ESI 

each, having normalized volumes 4 and 20 respectively. The corresponding 
facets of N{E2222) are given by inequalities of the form xqooo + a;oooi > 5. □ 

Example 22 (The largest tight span). There is a unique coarsest regular 
subdivision with 13 maximal cells. Its tight span is a 3-dimensional polytope 
with /-vector (13,24,13,1), as shown in Figure El The corresponding facets 
are given by inequalities like Xqooo + 2^1111 — a;ooii < 44. □ 

Example 23 (The largest tight span with non-pure dimension). Al- 
most half of the coarsest regular subdivisions have a pure dimensional tight 
span. The largest non-pure dimensional tight span has /-vector (11, 20, 11, 1). 
It contains one 3-dimensional face and two maximal 2-faces. A figure is de- 
picted in Figure ini The corresponding facets are given by inequalities like 

3^0000 + 3;ioio + s^iiii — a;iooo ~ 2^1001 ~ ^^loii < 47. |— | 
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All the 334 facet inequalities and all the data we have discussed in this paper 
are available on our website jbio . math . berkeley . edu/ 4cube/ 
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